Non-universal suppression of the excitation gap in chaotic Andreev billiards: 
Superconducting terminals as sensitive probes for scarred states 
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When a quantum- chaotic normal conductor is coupled to a superconductor, random-matrix theory 
predicts that a gap opens up in the excitation spectrum which is of universal size E^ MT « 0.3h/tD, 
where tD is the mean scattering time between Andreev reflections. We show that a scarred state of 
long lifetime ts ^> tD suppresses the excitation gap over a window AE w 2Ef MT which can be much 
larger than the narrow resonance width Ts = h/ts of the scar in the normal system. The minimal 
value of the excitation gap within this window is given by Ts/2 <C E^ MT . Hence the scarred state 
can be detected over a much larger energy range than it is the case when the superconducting 
terminal is replaced by a normal one. 

PACS numbers: 05.45.Mt, 03.65.Sq, 73.21.-b, 74.45,+c 
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Typical wavefunctions in classically chaotic quantum a ) 
systems [J Q display a structure-less speckle pattern 
which is well captured by Berry's random plain-wave 
model 3]. These quantum systems obey a strong degree ^ 
of universality which makes them amenable to a random- S 
matrix theory (RMT) description. Recent advances of 
semiclassical techniques 0, EI El link the random- 
matrix universality to the condition of well-developed 
wave chaos. Two sources of deviations from this con- 
dition have been identified. The first source is rooted in 
the quasi-deterministic quantum-to-classical correspon- 
dence of the dynamics for times shorter than the Ehren- 
fest time ts jj El 0] The second source is the 
scarring of wavefunctions along the trajectories of peri- 
odic orbits 0,EliniElE3 

The complementary nature of both sources of non- 
universality becomes apparent when the system is opened 
up, so that typical states acquire a width of order 
E^h — h/tp (the Thouless energy), where tjj is the clas- 
sical mean dwell time in the system |TJ. Quantum-to- 
classical correspondence then induces anomalously short- 
lived states with resonance width much larger than E^h 
0. Scarred states, on the other hand, can acquire life- 
times ts much larger than the typical lifetime t& when 
the openings are not visited by the periodic orbits which 
support the scar. The resulting narrow resonances of 
width Ts = h/ts have been observed experimentally and 
numerically in a large variety of physical systems, such 
as lateral quantum do ts llq . microwave cavities 17], res- 
onant tunnel diodes |l8|, micro-optical lasers 19], and 
surface waves [2(j. 

In general .scar red states are a rare fraction of all states 
of a system |$f |l0|, [ill |l2|, |l3| and their influence can be 
ignored for almost all energies, with the exception of the 
resonance intervals of size Ts- The purposes of this pa- 
per is to demonstrate that scarred states can be detected 
over much larger energy ranges when the normal open- 
ing is replaced by a superconducting terminal so that the 
system forms a so-called Andreev billiard 21]. Chaotic dow AE w 2E, 
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FIG. 1: Left panel: Excitation spectrum as function of the 
Fermi energy Ef, obtained by numerical computations for the 
Andreev billiard shown in the right panel. The lowest exci- 
tation energy (triangles) determines the excitation gap. The 
second, third, forth and fifth excitation energy are shown as 
filled circles. The energy of the scarred eigenstate is high- 
lighted by open squares. The dashed line shows the random- 
matrix prediction E^ MT for the gap. Other lines are guides 
for the eye. Right panel: The Andreev billiard is a segment 
of a Sinai billiard with W = 0.21 a, h = 0.58 a, and R = 0.8 a. 
The pair potential in the superconductor is Ao = 0.015 Ef- 
In the energy range of the left panel, there are N — 25 open 
channels in the superconducting lead. 

Andreev billiards have attracted much attention over the 
past years [22j because they display a gap E g in the ex- 
citation spectrum which is induced by the superconduc- 
tor due to the dynamical process of Andreev reflection 
(retro-reflections at the superconducting interface which 
convert electrons into holes or vice versa, at the expense 
of Cooper pairs which are absorbed or emitted from the 
superconductor). Random-matrix theory predicts that 
this gap is of size E^ MT « 0.3 E Th , where the Thouless 
energy refers to the system with the normal opening (the 
associated dwell time tn is identical to the mean time 
between successive Andreev reflections in the Andreev 
billiard) miH. 

We show that a long-lived scarred state with lifetime 
ts ^ try reduces the excitation gap over an energy win- 
RiYiT ^ y s% The scar hence determines 
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FIG. 2: The square moduli \u\ 2 of the electron and \v\ 2 
of the hole components of the wavefunction for E f /Et\i — 
707.434 < ^F,min (a), for E F = E F , min (b) and for E F /E T h = 
707.695 > E F , min (c). 

fundamental system properties over a much larger range 
of energies than in the open normal system. The smallest 
value of the excitation gap within the window is given by 
Eg, min ~ rs/2. The minimal size of the gap is there- 
fore much smaller than the RMT prediction £^ MT . In 
particular, the reduction of the excitation gap can be far 
stronger than the reduction due to both the quantum- 
to-classical correspondence for times smaller than the 
Ehrenfest time, observed in Refs. j^HH? and the meso- 
scopic fluctuations, studied in |2f| . 

We first present numerical evidence for these results 
and then develop a theoretical description which relates 
the excitation energy of scarred states in the Andreev 
billiard to properties of the scar in the normal billiard. 

The excitation energies e of an Andreev billiard are 
given by the eigenvalues of the Bogoliubov-de Gennes 
(BdG) equation 
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V 2 — Ep is the single-particle Hamilto- 



nian with the chemical potential in the superconductor 
set equal to the Fermi energy Ep in the normal part of the 
system. The wavefunction is composed of an electronic 
part u{y) and a hole part v(r) which are coupled via the 
superconducting pair potential A(r) with value Ao in the 
bulk of the superconductor. Typical Andreev billiards 
satisfy the separation of energy scales E^h < Ao <C ^f. 
One can then assume that A(r) = Ao is constant in the 
superconducting region and vanishes in the normal part 
of the system [22j . Solutions of Eq. JU come in pairs of 
excitation energies ±e, where one solution is of the form 
and the other solution is of the form ) . 
In order to convince ourselves of the visibility of scars 
in the excitation spectrum of an Andreev billiard we 



solved the Bogoliubov-de Gennes equation numerically 
for the chaotic billiard shown in the right panel of Fig. 
[251 127^ . The normal part is a de-symmetrized segment 
of a Sinai billiard. The superconducting terminal is at- 
tached to one of the straight walls. The left panel of 
Fig. H shows the positive excitation energies as a func- 
tion of the Fermi energy Ef, scaled in units of the Thou- 
less energy E^h (the dependence of the Thouless energy 
on the Fermi energy is taken into account in this rescal- 
ing). The second, third etc. excitation energies show the 
characteristic level dynamics of chaotic systems when one 
parameter is changing, and avoided crossings can also be 
observed. Random matrix theory predicts that the ex- 
citation gap (set by the the lowest positive excitation 
energy) fluctuates around the value E^ MT = 0.3£^xh 
(dashed line) [H 23]. However, there is an energy in- 
terval of width AE/Et\i ~ 0.26 in which the excitation 
gap drops distinctively below this value. The smallest 
gap -E^min = 0.0654£?Th is attained at Fermi energy 
^F,min = 707.567£ Th . 

Inspection of the wavefunction of the lowest eigenstate 
in this energy interval reveals an interesting feature: for 
Ep smaller than £?F,min [see Fig. Eta)], the electronic 
component u is nonuniform and shows clear signatures 
of scarring along a family of marginally stable periodic 
orbits which bounce between the straight segments of 
the billiard. Calculating the probability P e = j \u\ 2 dr 
(Ph = J \v\ 2 dr) of finding the quasiparticle in the elec- 
tron (hole) state in the normal part of the system, one 
finds P e = 0.915 and Ph = 0.078, i.e. the quasiparticle is 
dominantly electron-like (the remaining weight is in the 
superconductor). The next panel [Fig. Eft))] corresponds 
to the eigenstate at the minimum of the excitation gap 
in Fig. nja), showing that in this case both u and v are 
nonuniform, and the probabilities are P e = 0.520 and 
Ph = 0.474. In Fig. Etc) we show an eigenstate at a 
Fermi energy which is larger than £V,min- Here the hole 
component is scarred by the same family of periodic or- 
bits and the quasiparticle becomes dominantly hole-like 
with P e = 0.064 and P h = 0.929. 

These observations suggest that the suppression of the 
excitation gap and the scarring of the wavefunction stems 
from a long-lived scarred quasibound state of the normal 
open billiard. When the superconducting terminal is re- 
placed by a normal lead, a scarred state can indeed be 
detected for energies very close to £?F,min ~ see the left 
panel of Fig. El which shows a scattering state computed 
at the energy E = £V,mm- For a slightly different en- 
ergy E = 707.657£^Th, however, this scar is already very 
much diminished, as is shown in the right panel of Fig. 
El At this energy, the gap in the excitation spectrum of 
the Andreev billiard is still strongly suppressed. These 
numerical results demonstrate that a long-lived scarred 
state can be detected over an energy window much larger 
than its resonance width Ts in the normal billiard when 
the system is coupled to a superconducting terminal. 



3 




FIG. 3: Scattering wavefunction of the open normal dot at 
energies E — i^F,min (a) and E — 707.65722rh (b). The states 
are excited by an incoming wave in the third transverse mode 
of the lead. 

We now present a theory which explains the numerical 
observations. In particular, we relate the Fermi-energy 
dependence of the scarred excitation in the Andreev bil- 
liard to the properties of the scarred state in the open 
normal billiard. Such relations can be derived because 
the excitation spectrum of Andreev billiards can be cal- 
culated from the scattering matrix S(E) of the normal 
system using the quantization condition [2^ 

det[l + S*(E F - e)S(E F + e)) = 0, (2) 

which is valid in the limit e <C Ao <C E F . The scat- 
tering matrix is a unitary (and in absence of a magnetic 
field also symmetric) matrix of dimensions iVxiV, given 
by the number of open channels N in the opening. In 
Eq. (0 S(E F +e) [S*(E F -s)] describes the propagation 
of electrons [holes] and * denotes complex conjugation. 
Equation (0) embodies the condition of total construc- 
tive interference after two Andreev reflections, so that 
an electron is converted first to a hole and then back 
to an electron. (Each Andreev reflection contributes a 
phase factor of i.) 

Narrow resonances in the normal system are associ- 
ated to poles of the scattering matrix which are situated 
close to the real axis in the lower half of the complex 
energy plane. In what follows we will assume that there 
is only one such narrow resonance in the energy range 
of interest, namely that associated to the scarred state, 
located at a complex energy £ = Es — i In or- 

der to isolate the contribution of this narrow resonance 
to the scattering matrix we first diagonalize S(E) = 
0(E)A(E)0 T (E) where 0(E) is an N x N dimensional 
orthogonal matrix and A (22) is a diagonal matrix con- 
taining the eigenvalues Xi(E), i = 1 . . . N. In the vicin- 
ity of the narrow resonance the energy dependence of 
0(E) can be neglected 28] . Inserting the decomposition 
S(E) = OK(E)0 T into Eq. Q one then finds that the 
quantization condition can be written as 

Y[[l + \*(E F -£)\i(E F +e)}=0. (3) 

i 

One of the eigenvalues, \$(E), is associated to the narrow 
resonance, while the other eigenvalues are associated to 



the remaining non-resonant states in the system. The 
resonant eigenvalue can be approximated as Xs(E) = 
e 2i9 s {E) E ~£* [2^] where the energy dependence of 63(E) 
can again be neglected over the energy range of interest. 
Demanding that the term involving Xs in Eq. Q be zero, 

t E F -e s -Es-iTsl1E F + es-Es + iTsl* =Q 
E F - e s - E s + i T s /2 E F + e s - E s - i T s /2 

(4) 

one can find the excitation energies 

e s = ±y/(E F -E s ) 2 + T%/4 (5) 

associated to the scarred eigenstate of the Andreev bil- 
liard in terms of the energy and width of the scarred 
resonance in the open normal system. The excitation en- 
ergies come in pair as required by the particle-hole sym- 
metry of the BdG equation JIJ. 

For Fermi energies far less than the energy Es of the 
scar in the normal system, the positive excitation en- 
ergy is given by es ~ Es — E F , which is then far larger 
than the RMT gap and hence corresponds to a highly 
excited state. As E F is increased, the excitation drifts 
down through the spectrum, and assuming Ts <C E^ MT 
the scarred state eventually becomes the lowest excited 
state at an energy E F « Es - E^ MT . This defines the 
beginning of the energy interval in which the scar can be 
detected by measuring the excitation gap of the Andreev 
billiard. The minimal value E g ^ n = mm FF \es\ = 
of the excitation gap is obtained for 22p,min = Es- When 
the Fermi energy is further increased the excitation en- 
ergy increases as well, and for values E F > Es + E^ MT 
the scarred state becomes again a higher excited state, 
with energy es ~ E F — Es. This explains the qualitative 
features observed in Fig.^a). 

The microscopic parameters Es and Ts/2 entering Eq. 
(0 can be determined directly from properties of the 
open normal billiard. The left panel of Fig. 0] shows the 
phase Q(E) = Imlndet S(E) of the scattering matrix 
over the energy range in which the scar determines the ex- 
citation gap. The energy dependence due to the resonant 
eigenvalue Xs is given by ImlnA^ ~ 2 arctan 2 (e-e s ) 7 
while the non-resonant states in the system contribute 
a linear background CE + 60 with C « 2tt/5n, where 
Sn is the mean level spacing of the normal system. A 
fit of the curve 2 arctan 2 (^ r ?E s ) + CE + Qq to the nu- 
merical data delivers the values Es = 707.56722xh and 
r^/2 = 0.074522xh. The width is in good agreement 
with the smallest value 2^, m in of the excitation gap given 
above, and the energy of the scar coincides exactly with 
the Fermi energy 22p,min at which this minimal value is 
attained. The right panel of Fig. 0] shows that substitut- 
ing the obtained parameters Es and Ts into Eq. ©, a 
quantitative agreement can be found between the theoret- 
ical prediction and the computed values of the excitation 
energies of the scarred states. 
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FIG. 4: Left panel: The phase Q(E) = Im In det [£(£)] of 
the scattering matrix of the open normal billiard [Fig. Q^b) 
with the superconducting terminal replaced by a normal lead] 
is plotted as a function of energy (circles). The solid line 
shows a fit of the function 2 arctan 2 (e-e s ) + C ^ + ®o to 
O(E). Right panel: Quantitative comparison of the excitation 
energies of the scarred Andreev bound state (squares) with 
Eq. J5J (solid line), using the parameters Es and Ys of the 
scar in the normal system. 

Finally, we turn to the observation that the scarred 
character of the state in the Andreev billiard shifts from 
the electronic component for Ep <C -E^min to the hole 
component for Ep ^> ^F,min, while both components 
show an equal amount of scarring for Ep = ^F,min (see 
Fig.[2J). The weight of the scar in each component can be 
defined by the squared overlaps P e ,s = \ J ^s u dr\ 2 and 
Ph,s — 1/ ^gvdr\ 2 with the scarred wavefunction ips 
in the normal billiard. To a good approximation, these 
weights should be proportional to the time spent as a 
particle or hole excitation, which can be quantified by 
the Wigner delay time r w (E) = fi§ Q. Within our 
previous approximations, we can write 

^)= {E _ E J +T y, +KC. (6) 

Ignoring the constant C which accounts for the back- 
ground of non-resonant states we find 

P e ,s (E F -e s - E s ) 2 + r|/4 _ e s - (E F - E s ) 



h,S 



(E F +e s - Es) 2 + T|/4 £s + (E F - E s ) ' 



(7) 

For Ep <C Es this delivers P ey s ^> Ph,s, hence scarring 
in the electronic part of the wavefunction, while the hole 
part is scarred for Ep ^> Es- For Ep = Es the weights in 
both components are predicted to be equal. This explains 
the observations in Fig. [2j 

In summary, we have demonstrated theoretically and 
numerically that long-lived scarred wavefunctions of 
quantum-chaotic system can be observed over large 
ranges of the Fermi energy when the system is connected 
to a superconducting terminal. The resulting suppression 
of the excitation gap should be accessible in experiments 
by tunnel spectroscopy The Fermi energy can be 

controlled by side gates [3(| . 
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